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We consider the motion of ballistic electrons within a superlattice miniband under the influence 
of an alternating electric fleld. We show that the interaction of electrons with the self-consistent 
electromagnetic field generated by the electron current may lead to the transition from regular to 
chaotic dynamics. We estimate the conditions for the experimental observation of this deterministic 
chaos and discuss the similarities of the superlattice system with the other condensed matter and 
quantum optical systems. 

I. INTRODUCTION 

Recent progress in the fabrication of semiconductor superlattices has made investigations of the physical conse- 
quences of ballistic electron transport in these systems an area of considerable and growing interest. Studies of the 
interaction between ballistic electrons in superlattices and the electromagnetic field have already established that 
nonlinear phenomena may appear even for small field strength [1]. Theoretical investigations in this area have exam- 
ined oscillations of electrons due to their interaction with constant [2] and alternating [3] electric fields, self-consistent 
electron-field oscillations [4,5], propagation of electromagnetic solitons through the superlattice [6], and other prob- 
lems. Reviews of the early theoretical investigations of nonlinear effects in superlattice - field interactions have already 
appeared [7,8]. Recently, oscillations of electrons in superlattices under the infiuence of constant electromagnetic field 
(Bloch oscillations) were demonstrated experimentally [9] . 

It is well known, that for most of nonlinear dynamical systems with more than 1.5 degrees of freedom, transitions 
from the regular to the chaotic motion can take place [10]. Among quantum mesoscopic systems, for example, the 
possibility that nonlinear effects can lead to instabilities and, in particular, to dynamical chaos has been widely studied 
in Josephson junctions under the presence of an rf-field in a variety of papers involving theoretical modeling, numerical 
simulations, and experiments (for a review see, e.g. [11]). Several articles have also studied the possible occurrence 
of chaotic dynamics in superlattice-field interactions [12-15]. In particular, chaotic motion of ballistic electrons in a 
2D superlattice in a constant magnetic field was studied within the classical approximation [12], and recently this 
approach has been used to explain the experiments on magnetotransport in antidot arrays [13]. Other studies suggest 
that the transition to dynamical chaos can occur for superlattice electrons under the infiuence of a constant magnetic 
field applied perpendicular to the axis of superlattice and for either electromagnetic waves of constant amplitude [14] 
or for electromagnetic solitons [15]. Chaotic Larmor oscillations of the superlattice electrons [14,15] appear, when the 
carriers are initially populated rather close to the top of the superlattice miniband. 

In this paper we discuss another possible manifestation of dynamical chaos in superlattices: namely, the transition 
to chaotic behavior for ballistic electrons moving through a superlattice and interacting with spontaneously generated 
self-consistent field and with the external alternating electric field. We consider the specifically the case in which this 
transition to chaos may occur for electrons with energies belonging to the bottom of superlattice miniband. 

II. BASIC EQUATIONS 

Consider the motion of electrons within the miniband of ID superlattice under the the infiuence of an external 
alternating electric field applied along the superlattice axis z. Neglecting the inter-miniband transitions, the electron 
dynamics for weak fields may be described semiclassically using the effective Hamiltonian approach [8,16] 

i7 = e(p) + C/(r,t), p = P + -A(r,i), (1) 
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where e(p) is the (mini)band energy, P is the canonical quasi-momcntum of the electron, e is the modulus of the 
electron charge, v is the electron velocity, and A(r, t) is the vector potential due to the electromagnetic field. The 
crystal momentum of the electron Tik is related to the canonical electron momentum P by 



hk = p = P + -A(r,t). (3) 
c 

The potential energy U (r, t) in (1) is connected with the alternating electric field applied along z direction 

Ee^tit) = Eocosiftt). (4) 
In (4) Q is the frequency of the alternating electric field. The motion of electrons produces a current density 

jz = ~eNv^, (5) 

where N is the number of carriers per unit volume. Consequently, a "self-consistent" field, A{r,t), is generated and 
can be described by the Maxwell equation 
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At the same time, according to (1) and (2) this field affects the electron motion. 

We assume a standard dispersion relation for electron motion within the miniband corresponding to the tight- 
binding approximation [7,16], so that 
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where to* is an effective electron mass, a is a superlattice period and A is a miniband width. 

Prom (l),(2),(4),(5)-(7) one can derive a set of equations describing the interaction of electrons with the self- 
consistent field and with the alternating electric field 
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Consider a sample with a characteristic size shorter than the wavelength of the field A^. For typical mesoscopic 
samples, this leads to fields in the microwave or far-IR domain. Then in a first approximation one can neglect the 
spatial derivatives in (8) and rewrite this set of equations as one equation describing a parametrically forced pendulum 
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The frequency lus characterizes Bloch oscillations when the self-consistent field is absent and is often called the 
Stark frequency [7,8]. As first predicted by Epshtein [4], in the absenee of the alternating external field {ojs = 0), a 
periodic energy transfer between the field and the electrons is possible with the characteristic frequency uje- We note 
that nonlinear oscillations of this type in different semiconductors with nonparabolic dispersion laws were previously 
considered by Vatova [17]. 
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III. SOME ANALOGIES WITH OTHER CONDENSED MATTER AND QUANTUM OPTICAL 

SYSTEMS AND THE TRANSITION TO CHAOS 



Before examining in detail the nonlinear dynamics implied by Eq. (9), let us consider some instructive analogies 
with other systems. Eq. (9) is exactly the same as the equation describing a Josephson junction subjected to an rf 
field in the non-dissipative limit where the feedback mechanism is formed by an external circuit. As noted early, it 
is well-known that in the perturbed Josephson junction various dynamical instabilities and even dynamical chaos can 
be observed [11]. Quite recently Dunlap et al. [18] suggested connecting a superlattice in series to a capacitor and 
subjecting the system to an rf field. The circuit produces the feedback mechanism, and as a result, the system [18] 
demonstrates nonlinear properties, including the possibility of converting ac frequency to dc voltage. 

In these two references, the feedback mechanism is extrinsic, in that it is formed by the external circuit. In contrast, 
in our present system, the feedback mechanism is intrinsic and is formed through the influence of the self-consistent 
field on the electrons. This situation is analogous to the quantum optical system which consists of an ensemble of 
2-level atoms interacting with a self-consistent electromagnetic field and with an external time-periodic field [19-21]. 
The strength of the external field is characterized by the Rabi frequency, which is an analog of the Stark frequency uig 
defined in (10). The self-consistent oscillations of the field have a characteristic frequency the so-called "cooperative 
frequency" - which is an analog of the Epshtein frequency oje (10), and is also proportional to square root of the 
number of particles N^/'^. The miniband width in the frequency units A/7i is equivalent to the 2-level transition 
frequency and, pursuing the analogy to the end, the value ea formally corresponds to the dipole moment of the 2-level 
transition. The nonlinear dynamics of the system "2-level atoms plus self-consistent field plus external field" in some 
reasonable approximations can also be described by the equation of forced pendulum [19-21]. As we shall see below, 
these previous studies of analogous systems analogies from the viewpoint of nonlinear dynamics and transition to 
chaos will provide useful insights into the a priori complicated superlattice electron dynamics. 

We now turn to a detailed description of the nonlinear dynamics governed by equation (9). Using an expansion in 
terms of Bessel functions, one can rewrite equation (9) as 

oo 

^ + cul J2 (-l)"'/n(G)sin($-Fnm) = 0, (11) 

n=—oo 

where G = tos/fl and J„(x) is the standard Bessel function. The case when the external field is absent {u>s = 0) 
was first considered by Epshtein in [4]. It was shown in [4] that if the self-consistent field is initially present, then 
the nonlinear energy exchange between the ballistic electrons and the field is governed by the pendulum equation. 
The case of high-frequency perturbation fl ^ uje has been also considered by Epshtein in [5]. In this case one can 
neglect all terms with n 7^ in the expansion (11). Thus the only effect of the high-frequency perturbation is the 
renormalization of the frequency u)e uje\J Jq{G). In our case, the natural initial conditions for the eq. (9) are those 
for which both the self-consistent field and its vector potential at t = are absent: £^sc(0) = ^sc(O) = 0. (Without 
loss of generality we assume ^sc(O) = const. = 0). Taking into account that E = (— l/c)A, we have ^(O) = $(0) = 0. 
These initial conditions correspond to the elliptic stable fixed point of a pendulum without perturbation (wg = 0). 
In the remainder of this paper, we shall consider only these initial conditions. In the absence of the external field, the 
self-consistent field is not generated and can't influence on the electron motion. But for LUg 7^ 0, the self-consistent 
field can be generated spontaneously, due to instabilities of the motion for some values of the perturbation frequency 
O. We first consider the case of G ~ 1 and retain in the Bessel expansion in (11) only the terms up to | n |< 2. Then, 
we have from (11) after linearization ($ <C 1) 
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$ = a;|Ji(G)sinm. (12) 



From (12) it is evident that at « w_e and O w {2loe)/1 (with I an integer), instabilities exist corresponding to 
linear and to parametric resonances. The strongest resonances occur at J7 « loe and at O « 2wb. Of course, the 
growth of the self-consistent field $ due to the linear instabilities will eventually saturate at an amplitude at which 
the previously neglected nonlinear terms become significant. 

Let us now return to equations (9), (11) and consider the case G 3> 1. This is the well-known problem of nonlinear 
resonance crossing, which was one of the starting points in the investigations of the Hamiltonian chaos [22] . In another 
physical context, analysis of the same mathematical problem using the method of nonlinear resonance overlap (the 
"Chirikov criterion") [22,23] showed [24] that at 
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K « » 1 (13) 

02VG ' 

the nonlinear dynamics becomes chaotic for the majority of initial conditions, and the maximal amplitude of the 
generated self-consistent field -B™"^ in the units of frequency is 

— I S™"^ 1 = 1 |~ on = ivs. (14) 

Therefore, the maximal amplitude of the chaotic field is of the same order as the amplitude of the external field 
I l^^""^ |~ i^o- In contrast, at ii" < 1 the nonlinear dynamics is regular, and the maximal amplitude of the generated 
self-consistent field is of the order of or less than the width of a single nonlinear resonance 

^ I 1 = 1 |~ oje. (15) 

As one can see, the basic equation (9) can be transformed to the equation of the periodically forced pendulum 

^ + u%sm^ = gsinCtt, (16) 

by the substitutions \1/ = $ — Gsinfit and g = u)s^- 

As we noted above, in the present physical context we should solve eq. (9) (or (16)) for the initial conditions 
$(0) = <i>(0) = 0. Previously, it has been shown ([20]), for the quantum optical analog of the model described by 

eq. (16), that at Q./uje ~ 1 dynamical chaos is possible even when at t = the self-consistent field is absent. This 
situation can be realized if the dimensionless perturbation parameter g j lo\ is larger than some critical value of order 
one. 

So, for slow external perturbation {^/uje ~ 1) the self-consistent field can be generated spontaneously and be 
chaotic. It should be noticed that in this case the amplitude of the self-consistent field is of the same order as or 
larger than the field amplitude for the regular motion. 

The dynamical behavior of the system (9) is illustrated in Figs. 1-4. The chaotic time-dependence of the sponta- 
neously generated self-consistent field is shown in Fig. la for Q/cue = 0.5. In contrast, for the same amplitude of 
the external field, but for Q/oje ~ 1, the dynamics of the self-consistent field is regular (Fig. lb). It is seen from 
comparison of Figs, la and lb that the amplitude of the chaotically generated self-consistent field is several times 
larger than the field amplitude of the regular motion. Fig. 2 demonstrates the regions of regular motion and of 
strong chaotic dynamics in the plane of dimensionless parameters VL/uje and ws/uje, and at fixed initial conditions: 
$(0) = $(0) = 0. Figs. 3 and 4 illustrate the modification of the nonlinear oscillations of the self-consistent field 
under the variation of the amplitudes of the external field but at fixed frequency ^I/loe- Fig. 3 shows the nonlinear 
dynamics at rather slow frequency, O/wb = 0.1. In this case, the chaotic dynamics exhibits intermittent behavior (see 
Fig. 3a). In the region of high-frequency external field (VI/loe ^ 1), the transition from regular to chaotic dynamics 
is rather sensitive to variations of the parameters (see Fig. 2). The chaotic dynamics in this case also reveals the 
character of intermittency, which is demonstrated in Fig. 4a. 



IV. CONCLUSION 



Our analysis establishes that the motion of ballistic electrons through a semiconductor superlattice can, when 
one takes account of the generation of a self-consistent field, demonstrate both linear and nonlinear instabilities and 
deterministic chaos. We have argued that this behavior is quite analogous to the quantum optical model [20] describing 
nonlinear dynamics of 2-level atoms interacting with the self-consistent and with the external time-periodic fields. 

It is important to examine the condition of validity of our approach. Prom [7,8], we conclude that this condition 
requires that all characteristic frequencies are much less then the miniband width. In our case, this condition takes 
the explicit form 

max{fl,u)E,i^s} A/h. 

For A ^ 10~2cV, ^ lO'^'^cm"^, a ^ 10^*'cm and eaEa ^ ImeV [3,7], the Stark frequency (ojs) is much less 
than A/h but comparable to the collective Epshtein frequency cue (10) and thus belongs to the terahertz region. In 
this case, the criterion for chaotic dynamics may be satisfied. So, one can see that this first general condition for 
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a transition to chaos in the case of balhstic electrons in superlattices are close to those needed for the observation 
of Bloch oscillations and related phenomena in superlattices [2,3,7-9,18]. For applications to realistic experimental 
systems, it is essential to understand and model the role of dissipation caused by the collisions of the ballistic electrons 
with impurities and phonons. In our present considerations, the neglect of these effects means that our model nonlinear 
dynamics actually represents Hamiltonian chaos, whereas in realistic systems, the dynamics would likely be dominated 
by dissipative effects and consequently attractors. Nonetheless, for parameters similar to those discussed above, we 
expect dissipative chaos and resulting strange attractors to emerge. A more detailed discussion of these dissipative 
effects, of analytic estimates of the boundary of the region of chaos, and of potential applications to experiments will 
be presented elsewhere [25]. 
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FIGURE CAPTIONS 



. 1 The dependence of the dimensionless amplitude of the self-consistent field £ = eaEsc/fiujE on the dimensionless 
time T — ujEt: (a) chaotic dynamics for ujs/uje ~ 0.8, and Q./uje = 0.5; (b) regular dynamics for ujs/uje = 0.8, 
and ^/uje = 2. For both cases, the initial conditions are $(0) — ^{0) — 0. 

. 2 The regions with strong chaotic and with regular dynamics in the space of the dimensionless parameters: the 
frequency of the external field (fi/cj^), and the amplitude of the external field (075/0; £;). The initial conditions 
are $(0) = $(0) = 0. 

. 3 The same as in Fig. 1 but for the parameters: (a) chaotic dynamics for ujs/'^e = 1-8, and H/uje — 0.1; (b) 
regular dynamics for lus/uje = 1-2, and fl/cuE = 0.1. 

. 4 The same as in Fig. 1 but for the parameters: (a) chaotic dynamics for lus/loe — 1-2, and ^I/lue — 1.5; (b) 
regular dynamics for lus/loe = 1-5, and ^I/lue — 1.5. 
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